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Abstract 

First-principles calculations based on density functional theory and the 
pseudopotential method have been used to investigate the influence of gradient 
corrections to the standard LDA technique on the equilibrium structure and 
energetics of rutile TiC>2 and SnC>2 perfect crystals and their (110) surfaces. 
We find that gradient corrections increase the calculated lattice parameters 
by roughly 3 %, as has been found for other types of material. Gradient 
corrections give only very minor changes to the equilibrium surface structure, 
but reduce the surface energies by about 30 %. 
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I. INTRODUCTION 



Most first-principles calculations on condensed matter are nowadays based on density 



has to be made to the exchange-correlation energy, and the vast majority of calculations 
employ the local density approximation (LDA). The basic assumption is that the exchange- 
correlation energy per electron at any point in the system is related to the electron density 
at that point in the same way as in a uniform electron gas, and density gradients are ignored. 
Some theoretical justification can be given for thisi~i, and in practice the LDA works well 
in a wide range of situations. However, its accuracy is not always satisfactory, particularly 
when energy differences associated with changes of bonding are needed, as in e.g. molecular 
dissociation or the adsorption of molecules at surfaces. Attempts to improve the situation 
by adding lowest-order corrections in powers of the density gradient are not successful, but 
important progress has been made recently by requiring that the dependence of the energy 
on the gradients satisfies certain physical requirements. This has led to various forms of 
generalized gradient corrections 



In the last few years, there has been a large amount of work on the influence of different 



there has been little work on the effect of GGC on the properties of partially ionic materials 
such as the oxides Ti02 and Sn02 treated here. The surface properties of materials like 
these are extremely important, because of their application as gas sensors and catalysts. We 
have recently reported a detailed study of the bulk and surface properties of SnOj^S, and 
we have initiated work on the interaction of molecules with the surfaces of both Ti02 and 
An understanding of GGC is of considerable importance in this general area. The 
goal of the present paper is to study the effect on the bulk and surface properties of Ti0 2 
and Sn0 2 of the two widely used GGC schemes due to Perdew and Wang0>0 and Becke 



functional theory (DFT)Era. This theory is formally exact, but in practi 




ice an approximation 
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and PerdewEH. 

II. TECHNIQUES 

The calculations are performed using the pseudopotential methodi'i, so that only the va- 
lence electrons are represented explicitly, the valence-core interaction being represented by 
non-local norm-conserving pseudopotentials, which are generated by first-principles calcula- 
tions on isolated atoms. Periodic boundary conditions are used, with the occupied electronic 
orbitals expanded in a plane-wave basis. The expansion includes all plane waves whose ki- 
netic energy h 2 k 2 /2m (k the wavevector, m the electronic mass) is less than a chosen cutoff 
energy E cut - The inclusions of gradient corrections within the pseudopotential plane- wave 
technique has recently been discussed in detail by White and BirdU, who show that a ro- 
bust and accurate calculation of the GGC exchange- correlation energy and potential can be 
achieved by summation on exactly the same real-space grid as would be used for the LDA. 
This technique has been used in the present work. 

The first-principles pseudopotentials in Kleinman-Bylander representation^ were gen- 
erated using the optimization scheme of Lin et aZ.il in order to reduce the required value 
of the plane- wave cutoff E cut . The pseudopotentials used in the GGC calculations were 
constructed consistently by including gradient corrections in the generation scheme. The Sn 
pseudopotential was generated using the 5s 2 5p 2 configuration for s- and p-wave components, 
and the 5s 1 5p a5 5<i a5 configuration for the d-wave. The core radii were equal to 2.1, 2.1 and 
2.5 a.u. for the s, p and d components respectively. The Ti pseudopotential was generated 
using the 4s L85 3<i 2 configuration for s and d waves and the 4s 1 4p°' 5 3<i - 5 configuration for the 
p wave, with core radii of 2.2, 1.5 and 2.4 a.u. for s, p and d waves respectively. The oxygen 
pseudopotential used in our LDA calculations was generated using the 2s 2 2p A configuration 
for the s and p waves and the 2s 2 2p 2 ' 5 3<i a5 configuration for the d wave, with a single core 
radius of 1.65 a.u. For the gradient-corrected oxygen pseudopotential, we have used the 
single configuration 2s 2 2p 3 ' 5 3d 0A5 and the same core radius. The use of a core radius of 
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1.65 a.u. means that there is an appreciable overlap of the oxygen and metal core spheres 
in the SnC>2 and Ti02 crystals, and in principle this could cause inaccuracies. However, 
direct comparisons of the present results with our earlier work on SnO^, which employed 
an oxygen pseudopotential with the smaller core radius of 1.25 a.u., show that any errors 
due to core overlap are very small. The calculations have been done using a plane wave 
cut-off E cut of 600 eV for SnC>2 and 1000 eV for Ti02- Our tests show that with these 
cut-offs the energy per unit cell is converged to within 0.2 eV, the convergence with respect 
to E cnt being not noticeably influenced by the inclusion of gradient corrections, even though 
the gradient corrected pseudopotentials are less smooth and regular than the LDA onesElS. 

The calculations were performed using the CETEP code!! (the parallel version of the 
serial CASTEP code!) running on the 64- node Intel iPSC/860 machine at Daresbury Labo- 
ratory. The code uses the band-by-band conjugate-gradient technique to minimize the total 
energy with respect to plane-wave coefficients. The LDA calculations were performed using 
the Ceperley-Alder (CA) exchange-correlation function^. 

For the ground states calculations Brillouin zone sampling is performed using the lowest 
order Monkhorst-Pack set of k-pointsi^, as in our earlier work on SnC>2& Electronic densi- 
ties of states (DOS) associated with the ground state were calculated using the tetrahedron 
methodic with fc-point sampling corresponding to 750 tetrahedra in the whole Brillouin 
zone. 



III. RESULTS AND DISCUSSION 

A. Perfect S11O2 and TiC*2 crystals 

The 6-atom rutile unit cell of Sn02 and Ti02 is characterized by the two lattice parame- 
ters a and c and the internal parameter u: the positions of the four oxygens are (±ii, ±ii, 0), 
(| ±u,| T«, !)■ The equilibrium structure has then been determined by relaxation with 
respect to the lattice parameters a and c and the internal parameter u. The equilibrium 
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values of these parameters both with and without gradient corrections are given in Table |. 

As usually happens, there is a tendency for the LDA to underestimate the lattice param- 
eter. This is especially noticeable for SnC>2, where there may also be an effect due to our 
treatment of the Ad shell as part of the core. The inclusion of gradient corrections tends to 
increase the lattice parameters, as has already been found for semiconducting and metallic 
systems@0il. The increase is 4 % or more for the Perdew-Wang GGC, and leads to results 
for a and c that are appreciably greater than experimental values. For the Becke-Perdew 
GGC, the increase is roughly 3 %. Both the c/a ratio and the u parameter are almost 
unaffected, and this suggests that the gradient corrections have the effect of an isotropic 
negative pressure, as pointed out by Seifert et alM. 

We have calculated the electronic DOS for the Sn02 perfect crystal using both LDA and 
the two GGC schemes, but the changes caused by GGC are very small and we do not show 
the results here. 



B. The Sn0 2 and Ti0 2 (110) surfaces 

Our calculations on the stoichiometric (110) surface of both materials have been done 
with the usual repeating slab geometry. The rutile structure can be regarded as consisting 
of (110) planes of atoms containing both metal (M) and oxygen (O) atoms, separated by 
planes containing oxygen alone, so that the sequence of planes is O - M2O2 - O - O - M2O2 
- O etc. The entire crystal can then be built up of symmetrical 3-plane O - M2O2 - O units. 
The slabs we use contain three of these units, and our repeating cell contains 18 atoms (6 M 
and 12 O). The perfect (110) surface consists of rows of bridging oxygens lying above a 
metal-oxygen layer. The vacuum separating the slabs has been taken wide enough to ensure 
that interactions between neighboring slabs are small. The width we use corresponds to two 
O - M2O2 - O units, and is such that planes of bridging oxygens on the surfaces facing each 
other across the vacuum are separated by about 6.8 A. 

The surface structure has been determined by relaxing the entire system to equilibrium, 
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and the calculations have been done with and without gradient corrections. As in our 
previous work on SnC>2 (110)11, and the work of Ramamoorthy et al. on Ti02 (110)§, we find 
displacements of the surface atoms of order 0.1 A, with 5-fold and 6-fold coordinated metal 
atoms (Mn and Mi) moving respectively into and out of the surface, in-plane oxygens (On) 
moving out and bridging oxygens (Oi) moving very little. The changes of the bond lengths 
between the surface atoms, including sub-bridging oxygens (Om) and the uppermost oxygens 
(Oiv) of the following O - M2O2 - O unit, for LDA and gradient corrected calculations are 
given in Table |1[ From these results, it is clear that gradient corrections have only a minor 
effect on the relaxed equilibrium structure. As we have already noted for the perfect crystal 
case, modifications of atomic structure with respect to LDA results are more pronounced in 
the PW scheme. 

We have calculated the surface formation energy in the standard way, by subtracting 
from the slab total energy (18 atoms) three times the energy of a 6- atom perfect crystal 
unit cell and dividing by the total surface area. We find that the relaxed surface energy of 
Sn0 2 (110) is 1.66 Jm" 2 in the LDA, 1.13 Jm" 2 in PW-GGC and 1.16 Jm" 2 for BP-GGC. 
The LDA result is close to the value of 1.50 Jm~ 2 reported earlieril. For Ti02 (HO), the 
values are 1.14, 0.82 and 0.84 Jm -2 respectively. Our LDA result for Ti02 (HO) is close to 
the value of 1.06 Jm -2 reported by Ramamoorthy et a/.0. Comparison of LDA and GGC 
results shows that gradient corrections have a substantial effect on the surface energies. 
For both GGC schemes, the surface energies are lowered by about 30% with respect to the 
LDA values, the difference between PW and BP being very small. This decrease of surface 
energy by GGC is consistent with the general tendency of gradient corrections to remove 
the systematic overestimation of electronic binding energy in the LDA. 

The electronic DOS of the Sn02 (HO) surface using the LDA and the two GGC schemes 
are compared in Fig. |I[ In order to separate out effects of electronic structure, all the calcu- 
lations are done at the equilibrium lattice parameters and the relaxed positions produced by 
the BP scheme. Overall, the differences between the three sets of results are small. However, 
there are significant differences at the top of the valence band and at the top of the 0(2s) 
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band. As we found in our previous worko there are peaks at the top of both bands due to 
surface states, these states being concentrated on the bridging oxygens. The effect of GGC 
is to reduce the intensity of the peak at the top of 0(2s) band. The effect on the intensity 
of the valence band peak is less systematic, since BP increases it but PW decreases it. The 
reason why these effects are interesting is that there appears to be no experimental evidence 
for the surface-state peak at the top of the valence band, so that the LDA predictions seems 
not to be consistent with experiment. The present results suggest the possibility that this 
inconsistency may be due to inaccurate treatament of exchange and correlation. 

IV. CONCLUSIONS 

Our calculations show that gradient corrections increase the lattice parameters of TiC>2 
and SnO"2 by ~ 4 % for the Perdew-Wang scheme and ~ 3 % for the Becke-Perdew scheme. 
These effects are similar to those reported previously for metals and semiconductors. For 
the surfaces we examined, gradient corrections have very little effect on the relaxed surface 
structure, but the surface energies are substantially reduced - by ~ 30 % in both the Perdew- 
Wang and the Becke-Perdew schemes. The effects of gradient corrections on the electronic 
DOS of Sn0 2 (HO) surface are very small, except at the top of the 0(2s) and 0(2p) bands. 
The changes we find at the top of the 0(2p) band may be relevant to apparent inconsistences 
between calculated and experimental results for the surface DOS in this region. 
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FIGURES 

FIG. 1. The electronic DOS of the Sn0 2 (110) surface for LDA (solid line), PW-GGC (dotted 
line) and BP-GGC (dashed line). The calculations are made at the equilibrium lattice parameters 
of the BP scheme. For presentation purposes, we have broadened the calculated DOS by Gaussians 
of width 0.5 eV. 
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TABLES 

TABLE I. Comparison of theoretical and experimental values of lattice parameters a and c 
and the internal coordinate u of Sn02 and Ti02- The theoretical values are calculated using the 
Ceper ley- Alder form of LDA (CA), and the Perdew-Wang (PW) and Becke-Perdew (BP) forms of 
GGC. Experimental values are from Ref. 40. 







a 


(A) 


c 


(A) 


c/a 


u 


SnC-2 


CA 


4.645 


(-1.9%) 


3.060 


(-4.0%) 


0.659 


0.307 




PW 


4.868 


(2.8%) 


3.183 


( 0.0%) 


0.654 


0.307 




BP 


4.809 


(1.5%) 


3.159 


(-0.8%) 


0.657 


0.307 




expt. 


4.737 




3.186 




0.673 


0.307 


TiC-2 


CA 


4.625 


(0.7%) 


2.911 


(-1.6%) 


0.629 


0.305 




PW 


4.781 


(4.1%) 


3.072 


(3.9%) 


0.643 


0.305 




BP 


4.747 


(3.3%) 


3.039 


(2.7%) 


0.640 


0.305 




expt. 


4.594 




2.958 




0.644 


0.305 


TABLE II. Comparison of calculated bond length modifications on Sn02 


(110) and TiC-2 


(110) 


with respect to 


the bulk values for LDA 


(CA), and 


two alternative 


forms of GGC ex- 


change-correlation. 


















SnC-2 






Ti0 2 








CA 


PW 


BP 


CA 


PW 


BP 


Oiv - Mii 




-4.5% 


-4.9% 


-4.2% 


-5.6% 


-5.5% 


-5.6% 


Oi - Mi 




-3.7% 


-4.0% 


-3.8% 


-4.9% 


-5.5% 


-5.5% 


0„ - M n 




-1.2% 


-1.2% 


-1.2% 


-0.9% 


-1.2% 


-1.2% 


Oin - Mi 




4.3% 


4.3% 


4.8% 


4.6% 


4.7% 


4.5% 


On - Mi 




2.7% 


2.9% 


2.9% 


2.3% 


2.8% 


2.8% 
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